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Abstract By analyzing Chandra X-ray data of a sample of 21 galaxy 
groups and 19 galaxy clusters, we find that in 31 sample systems there ex- 
ists a significant central (R^10h~l kpc) gas entropy excess (AK ), which 
corresponds to - 0.1 - 0.5 keV per gas particle, beyond the power-law 
model that best fits the radial entropy profile of outer regions. We also find 
a distinct correlation between the central entropy excess A^ and ^-band 
luminosity L K of the central dominating galaxies (CDGs), which is scaled 
as AA'o °c L l K 6±0A , where L K is tightly associated with the mass of the su- 
permassive black hole hosted in the CDG. In fact, if an effective mass-to- 
energy conversion-efficiency of 0.02 is assumed for the accretion process, 
the cumulative AGN feedback E^^ k - nM BB c 2 yields an extra heating 
of - 0.5 - 17.0 keV per particle, which is sufficient to explain the central 
entropy excess. In most cases the AGN contribution can compensate the 
radiative loss of the X-ray gas within the cooling radius (- 0.002 - 2.2 
keV per particle), and apparently exceeds the energy required to deviate 
the scaling relations from the self-similar predictions (- 0.2 - 1.0 keV 
per particle). In contrast to the AGN feedback, the extra heating provided 
by supernova explosions accounts for 0.01 - 0.08 keV per particle in 
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groups and is almost negligible in clusters. Therefore, the observed corre- 
lation between AK and L K can be considered as a direct evidence for AGN 
feedback in galaxy groups and clusters. 

Key words: galaxies: active — galaxies: clusters: general — X-rays: 
galaxies: clusters — (galaxies:) intergalactic medium 

1 INTRODUCTION 

Over the past few decades many observational and theoretical efforts have been de- 
voted to the studies of galaxy groups and clusters. As of today, however, some fun- 
damental astrophysical processes that determine the basic properties of these celestial 
systems are still poorly understood. For example, the scaling relations between X- 
ray luminosity (L x ), gas temperature (T x ), gas entropy (K), and total gravitating mass 
(M total ), which are predicted by the self-similar gravitational collapse scenario (e.g., 
Kaiser 1986; Navarro, Frenk & White 1995), are challenged by observed deviations 
in such a distinct way that non-gravitational heating sources are invoked to dominate 
the gas heating process in the central ^100 kpc (e.g., Ponman et al. 2003; Donahue 
et al. 2006; Sun et al. 2009). Currently, most of the efforts have been focused on the 
AGN heating of the inter-galactic medium (IGM) (see McNamara & Nulsen 2007 for 
a review), since it is estimated that powerful AGN outbursts may repeat per 10 6 - 10 s 
yr and release 10 58 - 10 62 ergs per outburst into the environment, and this amount of 
energy is sufficient to balance gas cooling and heating on group and cluster scales (e.g., 
Rafferty et al. 2006, 2008; Birzan et al. 2009). 

However, there also exist problems in the AGN feedback scenario. For example, 
the absence of X-ray cavities, a natural product of the AGN activity, has been reported 
in about 40% cool core systems (e.g., Cavagnolo et al. 2008; Birzan et al. 2009), and 
only about 10% of quasars are found to host powerful radio jets (e.g., White et al. 
2007), which indicate that cavity- and jet-related feedbacks might not be generic. In 
NGC 4051, the observed mass and energy outflow rates due to the AGN activity is 4-5 
orders of magnitude below those required for efficient feedback (Mathur et al. 2009). 
Also, Jetha et al. (2007) reported that no significant difference of gas entropy profiles 
between radio-loud and radio-quiet galaxy groups is found (see also Sun et al. 2009). 

If AGN feedback does dominate the gas heating history in central regions of galaxy 
groups and clusters, it should be proportional to the central gas entropy excess AK, i.e., 
^feedback 00 ^ (Voit & Donahue 2005), where AK is measured beyond a power-law 
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model that best describes the gas entropy distribution in the intermediate and outer 
regions. On the other hand, by studying AGN cavities embedded in the X-ray halos 
of galaxy groups and clusters, Allen et al. (2006) found a close relation between the 
AGN feedback energy E^S. and the Bondi accretion power Psondi = r\Mc 2 oc Mg H 
(?7 is the mass-to-energy conversion efficiency, M is the accretion rate, and M BH is the 
black hole mass; Bondi 1952), which indicates that the AGN feedback should scale 
with the SMBH mass in the form of £j^ ack x Ml u . And we know that M Bll is related 
to the galaxy's bulge luminosity Lb u i ge and thus galaxy's ^f-band luminosity L K via 
M B h 00 ^buige 00 L K (Marconi & Hunt 2003; Batcheldor et al. 2007). Given above 
relations, we expect a tight correlation between the central gas entropy excess and the 
galaxy's ^-band luminosity, i.e., AK oc L 2 K , which has never been reported in literature. 

In order to examine whether this correlation holds or not, we analyze the Chandra 
archive data of a sample of 21 galaxy groups and 19 galaxy clusters to measure the 
central gas entropy excesses, and then compare them with the ^-band luminosities of 
central dominating galaxies (CDGs). In §2 we describe our sample, Chandra obser- 
vations, and data reduction. In §3 we measure gas density and temperature, and study 
the central gas entropy excess against L K . In §4 and §5 we discuss and summarize 
our results, respectively. Throughout the paper, we adopt the cosmological parame- 
ters H = 71 km s" 1 Mpc" 1 , Q b = 0.044, Q M = 0.27, and Q A = 0.73. Unless stated 
otherwise, the quoted errors stand for the 90% confidence limits. 

2 SAMPLE, OBSERVATION, AND DATA REDUCTION 
2.1 Sample selection 

In order to achieve our scientific goal, we need to select bright, nearby galaxy groups 
and clusters whose central galaxy can be well resolved with the Chandra Advanced 
CCD Imaging Spectrometer (ACIS). To investigate the interaction between the central 
AGN and the IGM, the selected systems are limited to have only one bright central- 
dominating galaxy and show no significant merger signatures. Our sample consists of 
20 brightest nearby (z < 0.1) galaxy groups and clusters selected from the flux-limited 
ASCA sample of Ikebe et al. (2002), and 18 galaxy groups from two ROSAT group 
samples constructed by Mulchaey et al. (2003) and Osmond & Ponman (2004), respec- 
tively. Besides, we add the group NGC 3402 (z = 0.0153), which satisfies the sample 
selection criteria, and the giant AGN cavity cluster MS 0735.6+7421 (z = 0.216) for 
comparison. All the selected systems lie above the Galactic latitudes of +15°, and are 
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located outside the fields of the Magellanic Clouds. Some basic properties of the sam- 
ple members are summarized in Table 1. 



Table 1 The sample 



Group/cluster 


Redshift 


CDG 


log(7.B/L B ) ° 


log(L sr /L K , ) t 


ObsID 


Chandra 
Date 


observation 

Raw (net) exposure (ks) 


NGC 0383 


0.0173 


NGC 0383 


10.86 


1 1.71+0.01 


2147 


2000-11-06 


47.2(44.4) 


NGC 0507 


0.0170 


NGC 0507 


11.02 


11.76+0.01 


2882 


2002-01-08 


43.7(43.4) 


NGC 0533 


0.0181 


NGC 0533 


11.07 


11.75+0.01 


2880 


2002-07-28 


37.6(37.4) 


NGC 0720 


0.0059 


NGC 0720 


10.49 


11.25+0.01 


7372 


2006-08-06 


49.4(48.8) 


NGC 0741 


0.0185 


NGC 0741 


11.19 


11.84+0.01 


2223 


2001-01-28 


30.4(29.6) 


NGC 1407 


0.0057 


NGC 1407 


10.77 


11.46+0.01 


791 


2000-08-16 


48.6(48.6) 


NGC 2563 


0.0159 


NGC 2563 


10.59 


11.42+0.01 


7925 


2007-09-18 


48.8(47.9) 


NGC 3402 


0.0153 


NGC 3402 


10.76 


11.42+0.01 


3243 


2002-11-05 


29.5(29.3) 


NGC 4125 


0.0050 


NGC 4125 


10.70 


11.28+0.01 


207 1 


2001-09-09 


64.2(63.6) 


NGC 4261 


0.0068 


NGC 4261 


10.72 


11.38+0.01 


9569 


2008-02-12 


101.0(99.3) 


NGC 4325 


0.0254 


NGC 4325 


10.66 


11.30+0.02 


3232 


2003-02-04 


30.1(29.9) 


NGC 5044 


0.0082 


NGC 5044 


10.85 


11.37+0.01 


9399 


2008-03-07 


82.7(82.7) 


NGC 5129 


0.0233 


NGC 5129 


11.07 


11.65+0.01 


6944 


2006-04-13 


21.0(21.0) 












7325 


2006-04-14 


26.7(25.8) 


NGC 5846 


0.0063 


NGC 5846 


10.70 


11.45+0.01 


7923 


2007-06-12 


90.4(89.8) 


NGC 6269 


0.0353 


NGC 6269 


11.38 


11.95+0.01 


4972 


2003-12-29 


41.4(39.6) 


NGC 7619 


0.0116 


NGC 7619 


10.98 


11.54+0.01 


3955 


2003-09-24 


37.5(28.8) 


HCG 42 


0.0128 


NGC 3091 


10.88 


11.60+0.01 


3215 


2002-03-26 


31.7(31.7) 


MKW 3S 


0.0450 


NGC 5920 


11.02 


11.57+0.02 


900 


2000-04-03 


59.1(57.2) 


JvlfV W 4 




INLrL. 4U / 3 




1 1 did A 1 




ZuUl- 1 1 -24 


1t~\ A( 1(1 Eft 


MKW4S 


0.0286 


NGC 4104 


11.20 


11.97+0.01 


6939 


2006-02-16 


37.1(35.8) 


UGC 12064 


0.0166 


UGC 12064 


10.48 


11.12+0.01 


■1057 


2003-09-18 


29.2(21.4) 


2A0335+096 


0.0349 


PGC 013424 


11.21 


11.83+0.02 


919 


2000-09-06 


21.4(19.7) 


Abcll 0085 


0.0556 


PGC 002501 


11.23 


12.05+0.02 


904 


2000-08-19 


40.5(38.4) 


Abell 0262 


0.0161 


NGC 0708 


10.60 


11.61+0.01 


7921 


2006-11-20 


111.9(110.7) 


Abell 0478 


0.0900 


PGC 014685 


11.32 


12.04+0.03 


1669 


2001-01-27 


42.4(42.4) 












6102 


2004-09-13 


10.4(9.9) 


Abcll 0496 


0.0328 


PGC 015524 


11.18 


11.87+0.02 


4976 


2004-07-22 


75.1(63.6) 


Abell 0780 


0.0538 


3C218 


11.38 


11.70+0.03 


4970 


2004-10-22 


100.7(98.8) 


Abell 1651 


0.0850 


PGC 088678 


11.00 


12.42+0.04 


4185 


2003-03-02 


9.7(9.6) 


Abell 1795 


0.0616 


PGC 049005 


11.01 


11.93+0.03 


493 


2000-03-21 


21.3(19.6) 


Abell 2029 


0.0767 


IC 1101 


11.33 


12.31+0.02 


4977 


2004-01-08 


77.9(77.7) 


Abell 2052 


0.0348 


3C 317 


10.87 


11.87+0.02 


5807 


2006-03-24 


127.0(126.9) 


Abell 2063 


0.0354 


PGC 054913 


10.79 


11.67+0.02 


6263 


2005-03-29 


16.8(16.8) 


Abell 2199 


0.0302 


NGC 6166 


11.38 


11.90+0.01 


497 


2000-05-13 


21.5(19.2) 


Abell 2589 


0.0416 


NGC 7647 


10.95 


11.85+0.02 


7190 


2006-06-11 


53.8(53.4) 


Abell 3 112 


0.0750 


PGC 012264 


11.35 


12.15+0.03 


2516 


2001-09-15 


17.5(16.5) 


Abell 3558 


0.0480 


PGC 047202 


11.39 


12.17+0.02 


1646 


2001-04-14 


14.4(14.2) 


Abell 3571 


0.0397 


PGC 048896 


11.55 


12.08+0.01 


4203 


2003-07-31 


34.0(31.6) 


Abell 4038 


0.0283 


IC 5358 


10.95 


11.68+0.02 


4992 


2004-06-28 


33.5(33.5) 


Abell 4059 


0.0460 


PGC 073000 


11.32 


12.00+0.02 


5785 


2005-01-26 


92.4(92.1) 


MS 0735.6+7421 


0.2160 


4C +74.13 


11.11 


11.85+0.07 


4197 


2003-11-30 


45.9(44.9) 



Notes: "fi-band luminosities of the CDGs shown as log^^/Lr^o): which are calculated using the data 
drown from http://leda.unfv-lyonl.fr (Paturel et al. 2003). *2MASS A"-band luminosities of the CDGs shown as 
log(L;;7LK.oX which are calculated using the data drown from http://www.ipac.caltech.edu/2mass 



2.2 Observation and data reduction 

All the Chandra observations (Table 1) were performed with the ACIS instrument op- 
erated at a focal plane temperature of -120 °C, and all the selected galaxy groups and 
clusters were positioned close to the nominal aim point on CCD 7 (ACIS-S) or on CCD 
3 (ACIS-I). We employed the Chandra data analysis package CIAO v4.1 and the cali- 
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bration CALDB v4.1.2 to process the datasets in the standard way, by starting with the 
Level-1 event files. We kept events with ASCA grades 0, 2, 3, 4, and 6, removed all the 
bad pixels and bad columns, and executed corrections for both the charge transfer inef- 
ficiency (CTI) and time dependent gain. We selected regions located as far away from 
the galaxy groups and clusters as possible to extract the 0.3-10.0 keV and 2.5-7.0 keV 
lightcurves for front-illuminated CCDs (CCD 0-3) and back-illuminated CCDs (CCD 
5 & 7), respectively, and then examined if there exist occasional background flares, the 
intervals contaminated by which were excluded in our study. In the spectral analysis 
that follows, we mask all the X-ray point sources detected in 0.3 - 8.0 keV with the 
CIAO tool celldetect, and apply the spectra extracted from Chandra blank-sky fields 
as the background. When source free regions are available in the observations, we also 
have attempted to use the spectra extracted from these regions as the background; we 
find that the results obtained with either background sets are consistent with each other 
within the 90% error limits. 

3 SPECTRAL ANALYSIS AND RESULTS 
3.1 Model fittings 

In order to measure the spatial distribution of specific gas entropy, which is defined as 
K = kTtie 2/3 , where T and n e are gas temperature and electron density, respectively, 
we divide each galaxy group or cluster into concentric annuli, and study the spectra 
extracted therein with XSPEC v 12.4.0 by applying a model that consists of an APEC 
component to represent the gas emission and a power-law component to represent the 
emission from unresolved point sources, both subject to a common absorption due to 
the neutral hydrogen. We set the metal abundance free to vary and fix the absorption 
at the Galactic value (Dickey & Lockman 1990), except that in 2A 0335+096 and 
Abell 478 the absorption is left free, because in these two galaxy clusters significant 
absorption excesses were reported in previous works (Mazzotta et al. 2003; Sanderson 
et al. 2005). We adopt the projct model embedded in the XSPEC package to correct 
the projection effect. By applying the F-test, we find that in the central annuli of the 
NGC 383, NGC 741, NGC 1407, NGC 4261, and NGC 5129 groups the power-law 
component is required at 90% confidence level, while in outer regions of these groups 
the power-law component is negligible. In other sample systems, the contribution of 
the power-law component is negligibly small. 

Using the best-fit deprojected spectral parameters, we calculate the 3-dimensional 
azimuthally-averaged radial entropy distribution K(R) for all the sample systems and 




Fig. 1 Deprojected azimuthally-averaged radial distributions of gas entropy 
K(R) of the galaxy groups (M 50 o < 10 14 M Q ; a) and galaxy clusters (M 50 o > 
10 14 M Q ; b) in the sample, along with the theoretical prediction K cc R lA with 
an arbitrary normalization (dashed lines), which is ascribed to shock heating 
that occurred during spherical collapses (Voit & Donahue 2005). 



plot them in Figure 1. Following Donahue et al. (2006), we fit the obtained radial 
entropy profiles with a three-parameter expression 

K(R) = AK + K m ( : f, (1) 

where AK represents the central gas entropy excess beyond the best-fit power-law 
model for larger radii, K wo is the normalization of the power-law component at 100/z^ 1 
kpc, a is the power-law index, and R is the 3-dimensional radius. All of the best-fit 
parameters are listed in Table 2. In all systems the gas entropy profile is well fitted 
with a power-law model in outer regions. In nine systems the power-law model can be 
extrapolate towards the center, resulting in a good acceptable fit. In other 31 systems, 
however, there exists a significant central (R < 10/z^ 1 kpc) entropy excess beyond 
the best-fit power-law model for outer regions. The obtained central entropy excesses 
vary in the range of - 100 keV cm 2 . This entropy excess is usually mentioned as 
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the central gas entropy plateau, and is suspected to be produced by heating sources 
like AGNs (e.g., Voit & Donahue 2005; McNamara & Nulsen 2007). Note that, in 2A 
0335+096, Abell 85, Abell 496, Abell 3558, Abell 3571, and Abell 4059 clusters, the 
evidence of minor merger has been reported in literature (Werner et al. 2006; Durret 
et al. 2005; Tanaka et al. 2006; Rossetti et al. 2007; Hudaverdi et al. 2005; Choi et 
al. 2004), which possibly result in mixing of high and low entropy gas (e.g., Ricker 
& Sarazin 2001; McCarthy et al. 2007). For these clusters, the average slope of the 
entropy profiles (1.13 + 0.17) shows no significant systematic bias from those of the 
sample (1.07 + 0.28), which indicate that minor merger likely does not appreciably 
affect the gas entropy distribution of the sample groups and clusters (Ghizzardi et al. 
2010). 

3.2 Correlation between central gas entropy excess and CDG's K-band 
luminosity 

In order to investigate the possible relation between the central gas entropy excess and 
the AGN heating, we apply the near-infrared .fif-band luminosity L K of the sample CDG 
by adopting the apparent K-baaA magnitudes of the CDGs from the Two Micron All 
Sky Survey (2MASS) archived and converting them into K-b&nd. luminosities (Table 
1). In the calculation, we have corrected Galactic extinction (Schlegel et al. 1998), and 
applied the correction for redshift (i.e., ^-correction) as k(z) = -61og(l + z) (Kochanek 
etal. 2001). 

In Figure 2, we show the central gas entropy excess AK versus the i^-band lumi- 
nosity L K for the 40 CDGs in our sample. It can be seen that AK Q apparently increases 
from to 100 keV cm 2 as L K increases from 1 x 10 11 to 3 x 10 12 L K ,©, and AK Q shows 
a roughly uniform scatter of - 0.5 in logarithmic scale. We find that the correlation 
coefficient for AK and L K is 0.45, which indicates that the two parameters has a cor- 
relation at 99% confidence level. We also have calculated the Kolmogorov-Smirnov 
(K-S) statistic of AK Q and L K against a proper constant model, and find that the cor- 
relation has a probability of more than 90% if the null hypothesis is true. Using the 
bisector of ordinary least-squares regression (Isobe et al. 1990), which is suitable to fit 
data with large scatters, we fit the log AK -\og L K relation with a liner model of 

log(A£ ) = A + B [log(L*/L K>0 ) - 10.9], (2) 



1 See 



http://www.ipac.caltech.edu/2mass 
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and obtain A = -0.8 ± 0.3 and B = 1.6 + 0.4 (Fig. 2), where the errors are determined 
by performing the Jackknife simulation that has been repeated 100 times (Feigelson & 
Babu 1992). 



+ 




jj I I J I U I J I II I ■ I I I I I L 



11 11.5 12 12.5 

log (L K /L K G ) 



Fig. 2 Central gas entropy excess A^ v s ^-band luminosity L K for the 40 
CDGs of our sample (see also Table 1 and 2). The solid line shows the best-fit 
model log(A^o) = -0.8 ± 0.3 + (1.6 + 0.4)[log(L r /L K , Q ) - 10.9] (Eq. 2 and 
§3.2), which is determined by the bisector of ordinary least-squares regression 
(Isobe etal. 1990). 



Table 2 Deprojected spectral analysis using an absorbed 
APEC+POWERLAW model 



Group/cluster 


N H a 


aV 7 


a L 


K 100 " 


Tx' 


M 5W > 


«»ol s 




E AGN i 

Icodbuck 


E J'JV ., 
feedback 






(10 20 cm~ 2 ) 


(keV cm" ) 




(keV cm 2 ) 


(keV) 


10 14 M© 


t}r} koc) 
71 


(ergs s~') 


(keV/gys particle) 


(keV/nas particle) 


(keV/gas particle) 


NGC 0383 


5.39 


55+ 7 - 47 
u - -0.55 


0.75 + 0.05 


389 + 47 


2 14+'* 
z - -0.60 


0.72 


1 1 


41.06 + 0.05 


2.26 


0.0047 


0.002±0.0001 


NGC 0507 


5.24 


17.8 ± 6.8 


140 +0.31 
'' w -0.22 


179 ± 37 


1.55 ± 0.06 


0.42 


52 


42 40+0- 09 
,z - -0.06 


4.76 


0.0123 


0.069±0.0002 


NGC 0533 


3.10 


J - -1.69 
6 1 1+5.93 
°- -5.81 
446+380 


1.43 ±0.04 


588 ±50 


^ 72+0.65 


0.50 


25 


42 03+0- 09 
-o.n 


3.85 


0.0115 


0.024±0.0001 


NGC 0720 
NGC 074 1 
NGC 1407 


1 .54 
4.43 
5.42 


091 +0.19 
""-0.11 

1 07+ 010 
'"'-0.08 

1 36+°- 12 
1 ■ -0.09 


100 ± 44 
420 ± 70 
1440 + 340 


44+8:11! 

-0.21 

1 80+°- 23 

-0.21 
j 77+O.34 


0.05 
0.54 
0.30 


40 
24 
16 


40 77+0- 22 
+u./ /_ 021 

41 62+0- 06 

-0.09 

40 8 1 +°- 1 5 

-0.32 


13.29 
4.37 
3.14 


0.0385 
0.0139 
0.0102 


0.017+0.0001 
0.009±0.000 1 
0.003+0.0001 


NGC 2563 


4.23 


60+^ 
-0.60 


1.30 + 0.05 


900 + 100 


1 69+<™ 

'■"-111 


0.49 


17 


41.17 ± 0.07 


1 .66 


0.0039 


0.003+0.0001 


NGC 3402 


4.43 


239+1.38 

4 68+- : ^ 
-e -2.57 

7 97+0.97 
-0.99 
5.25+ 1 - 77 


0.88 + 0.06 


122 ± 11 


0.93 + 0.04 


0.18 


^ 1 


42.49 ± 0.03 


5.12 


0.0179 


0.218±0.0002 


NGC 4 1 25 
NGC 4261 
NGC 4325 


1.84 
1 .55 
2.24 


84+0.21 
-0.14 
1.51 + 0.06 

1.30 ±0.09 


83 + 33 
2480 + 360 
141 + 15 


40+0-20 
"-0.18 

1 43+0- 54 
'■ -0.21 

1 07+0- 03 
1A "-0.02 


0.05 
0.37 
0.23 


33 
12 
58 


40 55+0.23 
WJ -Q.36 

40 78+0- 07 
4U/B -0.09 

42 83+0- 03 
+z.57_ 004 


16.88 
2.04 
2.89 


0.0750 
0.0072 
0.0109 


0.012±0.0001 
0.002+0.0001 
0.370±0.0003 


NGC 5044 


4.93 


■i n«+o:8 


1.38 + 0.03 


208 ± 8 


1.30 + 0.02 


0. ^2 


54 


42.60 ± 0.01 


2. ^8 


0.01 17 


0.150±0.0001 


NGC 5129 
NGC 5846 


l.7fi 
4.26 


39.5 ± 10.8 

3 15+°- 55 
J - -0.54 


i ,o.+0.24 
1.18+_ 018 

1.33 + 0.07 


247 ± 55 
431 +63 


98+0- 08 
"■ -0.05 

1.01 + 0.04 


0.20 
0.21 


16 

22 


41 56+ 006 
4 '- -0.07 

41 50+0- 17 
-0.29 


8.41 

-1.69 


0.0329 
0.0131 


0.023+0.0001 
0.019+0.0001 


NGC 6269 


4.77 


0.0I+ 981 


0.63 + 0.08 


226 ± 26 


7 C7+0.28 
z " -0.36 


0.94 


3 1 


42.19 ± 0.05 


3 13 


0.01 13 


0.017±0.0001 


NGC 7619 


5.00 


93+8:9i 
-0.74 


1.02 + 0.02 


424 ± 2 1 


1.03 ± 0.03 


0.21 


16 


41 28+0- 10 
*' z -0.15 


5.76 


0.0245 


0.01 1 ±0.0001 


HCG 42 
MKW 3S 


4.8 1 
i.03 


11 q+10.6 
"■'-10.7 
25 9+ 4 - 2 


1 47+°- 42 
-0.23 
1.10 + 0.09 


5 1 8 + 270 
149 ± 9 


93+0.10 
"'"-0.12 
3.96 + 0.10 


0.18 
1.95 


21 
106 


41 48+0- 10 
-0.09 

44.02 ± 0.01 


8.06 
0.5 I 


0.0236 
0.0021 


0.021+0.0001 
0.508+0.0001 


MKW 4 


1.89 


l 45+ 2 - 02 
''* J -l.45 


0.81 ± 0.04 


223 ± 18 


2 00+0- 07 
-0.05 


0.64 


47 


42.76 ± 0.02 


3.50 


0.0100 


0.099+0.0001 


MKW 4S 


1 .69 


n 0l +4.45 


0.99 + 0.1 1 


492 + 137 


1 8I+0- 75 
101 -0.48 


0.54 


24 


42 OS+O" 3 
4Z - -0.04 


-1.69 


0.0139 


0.025+0.0001 


UGC 12064 


1 1 .8 


2, +17.0 
-2.20 


66+ 016 
-0.13 


244 + 46 


1.65+Of 


0.47 


27 


41 62+ 006 
-0.07 


0.80 


0.0032 


0.010+0.0001 


2A0335+096 


20.3(1.78) 


3.70 + 6.85 


1.43 ± 0.04 


141 ± 5 


4 77+0 1 :T<j 
,zz -0.13 


2.18 


131 


44 40+0- 03 
•™-Q.05 


0.88 


0.0029 


1.072+0.0012 


Abell 0085 


3.42 


14.1 ± 4.5 


1.09 + 0.08 


174+ 12 


6.32 + 0.17 


4.21 


120 


44 44+0.03 


0.75 


0.0014 


0.556+0.0008 


Abell 0262 


5.38 


82+ - 89 
u - -0.82 


1.09 + 0.03 


304 + 19 


2.39 + 0.05 


0.86 


47 


42 87+^ : ^ 
4Z ° -0.05 


1 .43 


0.0021 


0.09 1 +0.000 1 


Abell 0478 
Abell 0496 


28.5(1.51) 
4.57 


3 8I+ 2 - 32 
-2 27 

5.89 + j ; S 


1.06 + 0.05 
1.15 + 0.05 


150 + 7 
178 + 8 


7 50+0- 38 
'■ -0.36 

5.29 ± 0.09 


5.49 
3.17 


1 S2 
109 


45 17+0- 04 
-0.08 
44.12 ±0.04 


0.53 
0.63 


0.0013 
0.0018 


2.218+0.0039 
0.369+0.0003 


Abell 0780 
Abell 1651 


4.92 
1.81 


i2.6+;| 

53.0+3" 


1.09 ±0.04 
0.80 ±0.14 


107 ±4 
234 ± 38 


3 47+O.O5 
J - ""-0.06 
7 45+0.83 
'■* J -0.77 


1.57 
5.45 


155 
96 


44 47+0-02 
-0.03 
44.24 ± 0.02 


0.92 
1.42 


0.0062 
0.0006 


1.810+0.0012 
0.258+0.0001 


Abell 1795 


1.19 


20.7 ± 4.2 


1.32 + 0.10 


127 + 10 


6.30 ± 0.20 


■1.17 


148 


44 74+0.03 
-0.06 


0.55 


0.0009 


1.104+0.0015 


Abell 2029 


3.05 


15.3 + 2.6 


1.06 ±0.04 


169 ±6 


8.10 + 0.15 


6.28 


160 


45 02+ 03 
* JUZ -0.05 


0.93 


0.0011 


1.338±0.0015 


Abell 2052 


2.73 


g 37+0.48 

65+^ 
-0.65 


1.48 ±0.02 


210 + 8 


3.20 + 0.04 


1.38 


109 


43 95+0 04 
,;, -* J -0.06 


1 .65 


0.0022 


0.634+0.0009 


Abell 2063 


2»K 


0.37 + 0.12 


170+16 


3.70 + 0.15 


1.75 


90 


43.64 ±0.01 


0.75 


0.0014 


0.239±0.0001 


Abell 2199 


0.86 


4 68+ 2 - 42 


0.97 + 0.06 


181+9 


4.80 ±0.14 


2.70 


120 


44 22+0-07 
-0.14 


0.83 


0.0034 


0.557±0.0018 


Abell 2589 


4.14 


o oi+™- 7 

u - -0.01 


0.45 + 0.10 


208 ± 13 


3.64 ±0.10 


1.70 


75 


43.52 ±0.01 


1.21 


0.0021 


0.185±0.0001 


Abell 31 12 


2.61 


o 37+379 

-VI 


1.12 + 0.13 


155 ± 20 


5 56+ 30 
J ' J -0.25 


3.37 


138 


44.58 +0.01 


1.23 


0.0024 


0.984+0.0002 


Abell 3558 


3.89 


67 4+ 33 - 6 


1.21+0.46 


199 ±41 


7.69 ±0.41 


5.83 


86 


43.92 +0.01 


0.70 


0.0014 


0.114±0.0003 


Abell 3571 


3.71 


92 6+ 20 - 1 


0.94 + 0.12 


193 ± 17 


7.17 + 0.23 


5.21 


101 


44 26+0 ° 1 
* H Z -0.04 


0.63 


0.0024 


0.290+0.0003 


Abell 4038 


1.56 


50+^? 


0.73 + 0.10 


186+14 


3.24 + 0.11 


1.41 


78 


43.41 ± 0.01 


0.98 


0.0028 


0.180+0.0001 


Abell 4059 


1.10 


8.28+]*J 


0.96 ± 0.05 


182 + 6 


4 36+ 00 ' 
4 ■ -0.08 


2.29 


101 


43 94+0.08 
<f J - -0.22 


1.30 


0.0035 


0.349±0.0018 


MS 0735.6+7421 


3.49 




1.26 ±0.14 


112+17 


5 37+0.41 

JJ 4ni 


2.94 


125 


44.56 + 0.01 


0.65 


0.0014 


1.079+0.0003 



Notes: "Absorptions are fixed to Galactic values (Dickey & Lockman 1990), except for those of 2A 0335+096 and Abell 478, which are left free and shown in 
brackets (see §3.1 for details). ^Central gas entropy excesses, power-law indices, and normalizations of Eq. 1 (§3.1). c Mean gas temperatures are measured from the 
spectra extracted within 0.1 - 0.2;"5oo- For each system, the mean gas temperature and r$oQ follow the relation r$oQ = 391 X T^ 63 /E(z) kpc (Willis et al. 2005).^Total 
gravitating masses within ('500 (§4). ^Cooling radii, at which the cooling time equals the universe's age at the system's redshift. ''0.3 - 12.0 keV luminosities are measured 
within the cooling radius in logarithmic scale. 'Cumulative AGN feedback to heating the IGM (§4 and Fig. 3), which is assumed £j^2u j — lMm{C 2 and rj = 0.02. 
^Supernova feedback to heating the IGM (§4 and Fig. 3), which includes both type la and II supernova contributions. *X-ray Radiative loss since z = 2 (§4 and Fig. 3). 
which is calculated from the 0.3 - 12 keV luminosity within the cooling radius (R coo {). 
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4 DISCUSSION 

By analyzing the deprojected gas entropy profiles, we find that there exists a significant 
central gas entropy excess in 78% sample groups and clusters, which can be ascribed 
to the non-gravitational heating processes. The average central entropy excess ranges 
from 5.7 keV cm 2 to 19.3 keV cm 2 from groups to clusters, which corresponds to a 
gas energy excess of - 0. 1 - 0.2 and =* 0.3 - 0.5 keV per gas particle, respectively, 
when either isodensity or isobaric heating process is assumed (e.g., Lloyd-Davies et 
al. 2000). Because the observed central entropy excess is correlated to the CDG's 
^-band luminosity, two of the most possible heating sources are AGN activity and 
supernova explosions. Here we compare their contributions and determine which one 
dominates the gas heating progress and is thus responsible for the observed central 
entropy excess. In the calculation that follows, we adopt gravitating mass M 50Q and 
gas mass fraction / gas ,5oo that are determined at r 50 o (/500 is the radius within which 
the over-density is 500 with respect to the universe's critical density at each system's 
redshift) following M 500 = E(z)' l \0 UA0 (kT x /3.0 keV) 165 * 004 M Q (Sun et al. 2009) 
and / gas = 0.0708 (kT x /l.O keV) 022 (Sun et al. 2009; Pratt et al. 2009), respectively, 
where E(z) = >/(l + z) 3 0. M + Q. A , and T x is the mean gas temperature measured in 
0.1 -0.2r 500 (Table 2). 

AGN activity 

It has been estimated that one powerful AGN outburst can produce a stable entropy 
excess of - 10 - 30 keV cm 2 for ~ 10 8 yr (Voit & Donahue 2005). In order to obtain 
the cumulative contribution of AGN feedback that can be approximated by £^™ ack - 
r]M BB c 2 (e.g., Croton et al. 2006; Short & Thomas 2009), we calculate the masses of 
the central black holes hosted in the sample CDGs by applying the relation between 
the near-infrared J^-band luminosity L K of the host galaxy and the black hole mass 
M BH , i.e., log(M BH /M ) = 8.21 + 0.07 + (1.13 + 0.12) x [log(Z^/L K , ) - 10.9], which 
is presented in Marconi & Hunt 2003 and shows similar scatter ranges to the more 
usually used relation between the stellar velocity dispersion of the galaxy's bulge and 
M B h (Marconi & Hunt 2003; Batcheldor et al. 2007). We adpot a conversion efficiency 
of 77 = 0.02, because in the study of X-ray cavities in nine elliptical galaxies Allen et 
al. (2006) found that about 2% of the accreting energy is converted into the thermal 
energy of the surrounding gas. By assuming that each gas particle within r 500 of galaxy 
groups and clusters has been heated uniformly, we find that the cumulative contribution 
of AGN feedback varies between 0.5 - 17.0 keV per particle for our sample, which 
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Fig. 3 Estimated energy feedbacks to the IGM by AGNs (circle) and su- 
pernova explosions (cross), comparing to the average gas energy excess of 
- 0.1 - 0.2 and 0.3 - 0.5 keV per particle for galaxy groups and clusters 
(grey belts), respectively, the X-ray radiative loss since z = 2 (triangle), and 
the energy required to deviate scaling relations from self-similar predictions 
(dashed line; Bode et al. 2009; see §4 for details). 



decreases as the gravitating mass M 50 o increases (Table 2 and Fig. 3), and is apparently 
higher than the average gas energy excess (- 0.1 - 0.5 keV per particle). 
Supernova explosions 

Using CDGs' 5-band luminosities (L B ; Table 1) and observationally constrained 
explosion rates of type la supernova (SN la; Dahlen et al. 2004; Cappellaro et al. 2005) 
and adopting that the supernova explosions heat the surrounding gas with ~ 10 50 ergs 
per event (e.g., Thornton et al. 1998), we find that the contribution of SN la feedback 
to the IGM since z = 3 (i.e., about 10-11 Gyr ago) is of the order of - 10" 4 - \0 2 
keV per particle. On the other hand, we note that, although no type II supernova (SN 
II) explosions is detected in nearby elliptical galaxies (Cappellaro et al. 1999), the 
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feedback energy of SN II exploded in high redshift and in galaxy's starburst epoch 
possibly has been deposited and contributed to the gas heating (e.g., Bryan 2000; Wu 
et al. 2002). Both field supernova observations in z = 0.1 - 0.9 (Dahlen et al. 2004; 
Cappellaro et al. 2005) and galaxy star formation theories based on initial stellar mass 
functions (e.g., Kravtsov & Yepes 2000) indicate that the explosion ratio of SN II to 
SN la is - 2 - 3. This value agrees with what were obtained in the studies on IGM 
metallicity (e.g., Wang et al. 2005; De Grandi & Molendi 2009), which suggested that 
the explosion ratio of SN II to SN la is ^ 1-6, the larger scatter of which is probably 
due to the clusters' different merger history. By assuming that the feedback energy per 
SN II event is the same as that of SN la event (e.g., Woosley & Weaver 1986), and 
that the cumulative explosion ratio of SN II to SN la is 3, we estimate that the total 
contribution of supernova feedback to the IGM is 0.0006 - 0.08 keV per particle, 
which also decreases as the gravitating mass M 500 becomes larger (Table 2 and Fig. 3). 
The estimated supernova contribution is apparently lower than the average gas energy 
excess (- 0.1 - 0.5 keV per particle), and mostly is about 1-2 orders of magnitude 
lower than the average gas energy excess. For all the sample groups and clusters, it is 
about 2-3 orders of magnitude lower than the AGN feedback. 

In most case the estimated AGN contribution can compensate the radiative loss 
since z = 2 (i.e., about 9-10 Gyr ago) in galaxy clusters (- 0.1 - 2.2 keV per 
particle; Table 2 and Fig. 3), which is calculated from the 0.3 - 12 keV luminosity 
within the cooling radius (Table 2). For galaxy groups (< 10 14 M©), the AGN feedback 
energy is about 1-2 orders of magnitude higher than the X-ray radiative loss (^ 
0.002 - 0.4 keV per particle; Fig. 3). Moreover, this surplus energy fed by AGNs in 
galaxy groups is expected to re-distribute the IGM gas, especially in the central regions, 
and thus break the self- similarity between galaxy groups and clusters. Assuming a 
uniform baryon (mainly includes gas and stellar components) mass fraction for galaxy 
groups and clusters, the energy required to deviate the scaling relations from self- 
similar predictions is - 0.2 - 1.0 keV per particle (Bode et al. 2009), which can be 
supplied by AGN not by supernova explosions. 



5 SUMMARY 

In 31 galaxy groups and clusters we find that there exists a significant central gas en- 
tropy excess, which scales with the J^-band luminosity of the CDG via A^ 00 L l K 6±0A . 
By comparing between the contributions of AGN activity and supernova explosions, 
we conclude that AGNs are responsible for the central entropy excesses. 
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